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\Q , The optical/near-infrared (OIR) region of the spectra of low-mass X-ray binaries 

■ appears to lie at the intersection of a variety of different emission processes. In this 

paper we present quasi-simultaneous OIR-X-ray observations of 33 XBs in an attempt 
to estimate the contributions of various emission processes in these sources, as a func- 
tion of X-ray state and luminosity. A global correlation is found between OIR and 
X-ray luminosity for low-mass black hole candidate XBs (BHXBs) in the hard X-ray 
state, of the form £oir K ^x 6 - This correlation holds over 8 orders of magnitude in 
Lx and includes data from BHXBs in quiescence and at large distances (LMC and 
M31). A similar correlation is found in low-mass neutron star XBs (NSXBs) in the 
hard state. For BHXBs in the soft state, all the near-infrared (NIR) and some of the 
optical emission is suppressed below the correlation, a behaviour indicative of the jet 
switching off/on in transition to/from the soft state. We compare these relations to 
theoretical models of a number of emission processes. We find that X-ray reprocessing 
in the disc and emission from the jets both predict a slope close to 0.6 for BHXBs, 
and both contribute to the OIR in BHXBs in the hard state, the jets producing ~ 90 
percent of the NIR emission at high luminosities. X-ray reprocessing dominates the 
OIR in NSXBs in the hard state, with possible contributions from the jets (only at 
high luminosity) and the viscously heated disc. We also show that the optically thick 
jet spectrum of BHXBs extends to near the if-band. OIR spectral energy distributions 
of 15 BHXBs help us to confirm these interpretations. We present a prediction of the 
Lom~Lx behaviour of a BHXB outburst that enters the soft state, where the peak 
Lorn m the hard state rise is greater than in the hard state decline (the well known 
hysteretical behaviour). In addition, it is possible to estimate the X-ray, OIR and radio 
luminosity and the mass accretion rate in the hard state quasi-simultaneously, from 
observations of just one of these wavebands, since they are all linked through correla- 
tions. Finally, we have discovered that the nature of the compact object, the mass of 
the companion and the distance/reddening can be constrained by quasi-simultaneous 
OIR and X-ray luminosities. 

Key words: accretion, accretion discs, black hole physics, ISM: jets and outflows, 
X-rays: binaries 



1 INTRODUCTION 
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X-ray binaries are binary systems in which a compact ob- 
ject - either a black hole or a neutron star - accretes 
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matter from a companion star. The origin of the emission 
from these sources is known in some wavebands and not 
so well established in others. Radio emission is produced 
by the synchrotron pro cess in collimated outflows (see e.g. 
iHiellming fc Hanlll995fl . whereas X-ray emission could orig- 
inate e.g. directly from the hot inner accretion disc, from 
a Comp tonising corona, from an advective flow, or from 
the j ets jPoutanen||l998|: iN^xavMi^Mah^devan^^^uataer^ 
199£ : iMcConnell et all l2000l iMarkoff. Falcke fc Fendeil 



200ll:lBrocksopp. Bandvopadhvav fc Fendedl2004ft~ 

The optical/near-infrared (optical/NIR; OIR) is per- 
haps the waveband for which the dominating emission pro- 
cesses in XBs are least known. OIR emission has been ex- 
tensively studied in outburst and quie scence (for reviews see 
Ivan Paradiis fc McClintocklll995 | and ICharles fc Coell2006t 



see also IChen. Shrader fc Livid ll997T) . The complex vari- 
ety of spectral, timing and luminosity properties of the 
OIR emission indicates that many processes may be con- 
tributing, each depending on a number of factors. In high 
mass X-ray binaries (HMXBs), the OIR light is largely 
dominated by the massive companion star in th e system 
Jvan den Heuvel fc Heiselll972l: iTreves et a?1ll980l) with oc- 
casional additional contributions, for example from the re- 
processing of X-rays. For low-mass neutron star X-ray bina- 
ries (NSXBs), there is strong evidence for a central X-ray 
source illuminating a disc that reprocesses th e light to OIR 
wavelengths (see e.g. iMcClintock et aljll979T) . 

The optical emission of a low-mass black hole candi- 
date X-ray binary (BHXB) is generally thought to arise 
in the outer accretion disc as the result of X -ray repro- 
cessing llCunnin gham[^97^J^tilek^^^2^2Sll: much like 
the NSXBs fe.g. lvan Paradiis fc McClintocklll995ft . Indeed, 
timing and spec tral analysis in man y cases has led to this 
conclusion (e.g. |Wa^ier_et_alJ |l99l|; |CaIlMari_et_alJ Il995t 
lO'Brien et alJl2002l: iHvnes et al.ll2002al ; lHvnedl2005l) . How- 
ever, reprocessed X-rays are often misleadingly assumed to 
dominate the OIR light in BHXBs. Some observations of 
BHXBs point towards alternative physical processes con- 
tribu ting (and sometim es dominating) the OIR emission 
(e.g. iHoman et alJl2005r) . 

Intrinsic thermal emission from the viscously heated 
outer accretion disc is expected to contribute signifi- 
cant light in the optical, through UV to X-ray wave- 
lengths JShakura fc Sunvaed Il973i: iFrank. King fc Raind 
2002). OIR behaviour has revealed this process to play a role 
in some BHXBs fe.g.lKuulkerdllll^lSoria Wu fc Johnston 



19991: IBrocksopp. Groot fc Wilmd l200ld : iBrocksopp et aT 
200ld : lHoman et al.l2005l) . Thermal emission from the com- 
panion star i s observed in some low - mass XBs in qu i- 
escence (e.g. lOkd Il977t iBailvnl Il992t lOrosz et all Il99 



iGreene. Bailvn fc Orosd 120011: iMikoIaiewska et alJ 12005ft . 

In the last decade evidence has been mounting for the 
flat optically thick spectrum of the jets to extend from 
the radio to the OIR regime <|Han fc Hjellmin 
Fenderl liool ICorbel fc Fenderl 120021: IMarkoff et al 
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1992; 
2003; 



Chatv et al.ll2003l: IBrocksopp et alJ 12004 iBuxton fc Bailvn 



2004 lHoma.n et al] 1 200.4 ~a~ view see i Fenderl 

2006). Behaviour that is not consistent with intrin- 
sic disc or reprocessed emission has in the past been 
attributed to e. g. magnetic loop reconnection in the 
outer disc (e.g. IZurita. Casares fc Shahbad 12003ft . emis- 
sion from a magnetically dominated compact corona (e.g. 



iMerloni. Pi Matteo fc Fabianl20o"o |) or emission from an ad- 
vective region fe.g. IShahbaz et al.ll2003ft . 

In BHXBs (both transient and persistent), proper- 
ties of the emission in all wavebands are often related 
to changes in the X-ray spectrum. The two main X-ray 
spectral states are the hard (or low/hard) state, which is 
characterised by a hard power-law spectrum and strong 
variability and the soft (or high/soft; thermal- dominant) 
state, where a thermal spectrum dominates with a pow er- 
law contribution (see IMcClintock fc Rcmillardl [200d for 
a review of X-ray states; see also ^ ^^^^^^^T l200ll : 
iFender. Belloni fc Galldl2004l: IHoman fc Bellonill200d) . The 

low luminosity 'quiescent' state is likely to be an extension to 
1 11 1 1 1 

the hard state !|NaraY anl|l996l:|Esin 1 McClintock fc Narayanl 

[ 19971: IMcClintock et aLlboodlFender Gallo fc Jonkerll2003l : 
IFender et al.ll2004t iGallo et al.ll200eh but currently this is 
not universally accepted. Here, we treat quiescence as an ex- 
tension to the hard state but we also show how our results 
differ when the quiescent data are removed. We hereafter 
class 'optical' and 'NIR' emission as that seen in the BVRI 
(~ 4400 - 7900A) and JHK (~ 1.25 - 2.22/xm) wavebands, 
respectively. 



1.1 Towards a Unified Model for the OIR 
Behaviour in BHXBs 

Power-law correlations between OIR and X-ray luminosities 
are naturally expected from a num ber of emission processes. 
Ivan Paradiis fc McClintock! |l994) showed that the optical 
luminosity of an X-ray reprocessing accretion disc varies as 
Lopt oc T 2 oc L x 5 a, where T is the temperature and a is the 
orbital separation of the system, and that this correlation 
has been observed in a selection of low-mass XBs. LoiR-ix 
correlations are also expected when the OIR originates in 
the viscously heated disc as both X-ray and OIR are linked 
through the mass accretion rate (see Section 3.2). 

In addition, OIR-X-ray correlations can be pre- 
dicted if the OIR emission originates in the jets. Mod- 
els of steady, compact jets demonstrate that the total 
jet p ower is related to the rad i o luminosity as L ra dio oc 
Lht llBlandford fc KonigH 119791: iFalcke fc Biermand Il996l: 
IMarkoff et alJl200lft . It was shown that the jet power is 
linearly proportional to the mass accretion rate in NSXBs 



and BHXBs in the hard state fFalcke fc Biermaml ^9' 
iMigliari fc Fendeil liool iKording. Fender fc Migliaril 1200 
and the X-ray luminosity scales as Lx oc m and Lx oc 
m 2 for rad iatively efficient and inefficient objects, respec- 
tively fe.g IShakura fc Sunvaevl Il973t iNaravan fc Yilll995l : 
IKording et al.ll200d) . The accretion in hard state BHXBs 
is found to be radiatively inefficient (the majority of the 
liberated gravitational potential is carried in the flow and 
not radiated locally), where jet-dominated states can ex- 
ist, whereas in NSXBs, the accretion is radiatively effi- 
cien t, and jet-dominat ed states are unlikely to exist (see 
also IFender et af]l2003l) . We therefore have: 



BHXBs: L ra dio <x L jct oc m 
NSXBs: L radio oc Lj c f oc m 14 



ocLx' 7 
ociV 4 



The correlation for BHXBs has been o bserved llCorbel et alJ 
| 2003t IGallo. Fender fc Pooled 12003ft and very recently, 
lMi^h^r^^^n^TT2^06l) have applied this technique to 
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NSXBs and found L ra dio oc L^ ; which is also consistent 
with the above NSXB model. If the optically thick jet 
spectrum is indeed flat from the radio regime to OIR, we 
can expect the following correlations: 

BHXBs: Loir oc L rad io oc L x - 7 

NSXBs: Loir oc L ra dio oc L x 4 

iHoman et alJ (I2005T) discovered a correlation between 
the quasi-simultaneous NIR (which was shown to originate 
in the jets) and X-ray fluxes for GX 339-4 in the hard state, 
with a slope Fnir oc f£ 53±0 - 02 ( 3-100 keV). To date, no 
other sources have been tested for jet OIR emission using 
OIR-X-ray correlations. 

It is now becoming clear that this profitable but sim- 
ple technique of analysing the dependence of OIR and X- 
ray luminosities over many orders of magnitude, may prove 
fruitful for the understanding of the emission mechanisms in- 
volved. The unification of jet-X-ray state activity is now un- 
derway; a steady jet exists in the hard state, which is acceler- 
ated as the X-ray spectrum softens, and is fin ally quenched 
as it passes the 'jet line' into the soft state dFender et alJ 
2004). A unification (if one exists) of the origins of OIR light 
from BHXBs in different spectral and luminosity states is 
desired to understand the behaviour of these systems. Fur- 
thermore, a measure of the level of OIR emission from jets 
may be used to constrain jet power estimates. 



2 METHODOLOGY & RESULTS 

For this work, we have collected OIR and X-ray data from 
a large number of BHXBs, NSXBs and HMXBs in order 
to find relations that may help determine the processes re- 
sponsible for the OIR light in these systems. We apply the 
technique of testing the dependency of OIR luminosity with 
X-ray luminosity for the three types of XB and between 
different X-ray states in BHXBs, and attempt to identify 
the dominant emission mechanisms for BHXBs at a given 
luminosity and X-ray state. 

A literature search for quasi-simultaneous (no more 
than ~ I day between observations; for sources with out- 
bursts <1 month in length we only use data with separa- 
tions of <0.1 times the outburst length) X-ray and OIR 
fluxes from BHXBs was conducted. Where possible, tabu- 
lated fluxes or magnitudes were noted. In some cases we 
obtained data directly from the authors. We also made use 
of the DEXTER applet provided by NASA ADS to ex- 
tract data from light curves where the data themselves were 
unattainable. For each source, the best estimates of its dis- 
tance, optical extinction Ay and HI absorption column TVh 
were sought. Table 1 lists the properties of each BHXB for 
which data were obtained. We used non-simultaneous OIR- 
X-ray luminosities only in quiescence for some sources, and 
for these we have included errors that encompass all ob- 
served values of the quiescent flux in one of the two wave- 
bands. 

The X-ray unabsorbed 2- 10 keV flux was calculated for 
all X-ray data . We made note of the X-ray state of each 
source on each observation, as defined by the analysis by 



authors in the literature. We were unable to apply strict 
definitions to the spectral states due to the differing nature 
(e.g. X-ray energy ranges and variability analysis) of each 
data set. We have therefore used the judgement of the au- 
thors to determine the spectral states; a method likely to be 
much more accurate than any conditions that could be im- 
posed by us. We assume a power-law with a spectral index 
a = —0.6 (photon index Y = 1.6) where F v oc v a , when the 
source is in the hard state, and a blackbody at a temperature 
of 1 keV for soft state data. These model s for the X-ray spec - 
trum are the same as those adopted bv lGallo et alJ J2003F) , 
and altering the values in the models to other reasonable 
approximations does not significantly change the estimated 
X-ray luminosities. For GRS 1915+105, we use data from the 
radio-bright plateau st ate, which is approximat ely analogous 
to the hard state (e.g. iFender fc Poolevll2000l) . The NASA 
tool Web-PIMMS was used to convert from instrument X- 
ray counts per unit time (e.g. day- averaged RXTE ASM 
counts s~ 1 ). lBrocksopp et alJ i2004h also provide a table of 
approximate instrument counts-flux conversion factors used 
in this work. 

For the OIR luminosities, data was collected from the 
optical B-band at 440 tim to the near-infrared L"-band 
at 2220 nra. OIR absorbed fluxes were de-reddened using 
the best-known value of the extinction Av to each source 
and the dependence of extinction with wavelength given by 
ICardelli. Clayton fc Mathisl l|l989h . For OIR data at fluxes 
low enough for the companion star to significantly contribute 
(i.e. quiescence for most low-mass XBs), the data was dis- 
carded if the fractional contribution of the companion had 
not been estimated in the OIR band of the data. The frac- 
tional contribution of the companion (column 7 of Table 
1) was subtracted from the low-flux OIR data in order to 
acquire the flux from all other emission processes in the 
system. This contribution is not well constrained in many 
cases due to the uncertain spec tral type of the compan- 
ion star (e.g. Haswell et al. 2 0021) . We have therefore prop- 
agated the errors associated with this into the errors of the 
OIR luminosities for all quiescent data. In sources for which 
the companion is c omparatively bright (e.g. GRO J1655-40; 
iHvnes et al.l ll998l. its contribution has been subtracted in 
outburst in addition to quiescence. 

The intrinsic (de-reddened) OIR and X-ray luminosities 
were then calculated given the best-known estimate of the 
distance to each source (Table 1). We adopted the approx- 
imation Loir ~ vF v ,oir to estimate the OIR luminosity 
(we are approximating the spectral range of each filter to 
the central wavelength of its waveband). The errors associ- 
ated with the luminosities are propagated from the errors 
quoted in the original data. Where no errors are quoted, we 
apply a conservative error of 30%. Errors associated with 
estimates of the distance, extinction Ay and NI absorption 
column JVh were sought (Table 1). Where these limits are 
not directly quoted in the reference, we used the most con- 
servative estimates implied from the text. These errors were 
not propagated into the errors associated with the luminosi- 
ties we derive because the resulting plots would be domi- 
nated by error bars, however we also show one error bar for 



This energy range was adopted to be consistent with that used 



for the radio-X-ray correlations of lGallo et"aH (2003; 2-11 keV) 
and lMigliari fc FendeJ (2005; 2-10 keV). 
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Figure 1. OPT and NIR refer to the optical (BVRI) and NIR (JHK) wavebands, respectively. As a rule, filled symbols represent NIR 
data and unfilled symbols represent optical data. The mean la errors associated with each value (which include errors on distance and 
reddening) are shown top left. Panel (a): Quasi-simultaneous X-ray versus OIR luminosities for BHXBs in the hard X-ray state. The 



best power-law fit is Lqir = 10 



13.1±0.6 r 0.61±0.02 



The black dashed lines represent the la uncertainties in the normalisation of the fit. 



Panel (b): As panel (a) but for BHXBs in the soft state (the BHXB hard state correlation is shown for comparison). 
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Low-mass NSXBs in the hard state 
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Figure 2. OPT and NIR refer to the optical (BVRI) and NIR (JHK) wavebands, respectively. As a rule, filled symbols represent NIR 
data and unfilled symbols represent optical data. The mean lcr errors associated with each value (which include errors on distance and 
reddening) are shown top left. Panel (a): Quasi-simultaneous X-ray versus OIR luminosities for NSXBs in the hard X-ray state. The 



best power-law fit (green dashed line) is -Loir = 10 



10.8±1.4 r 0.63±0.04 



The BHXB hard state correlation (red line with black dashed lcr 



uncertainties) is shown for comparison. Panel (b): As panel (a) but for HMXBs (black hole and neutron star XBs). 
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Table 1. Properties and data collected for the 16 BHXBs. 



Source 


Distance 


Period 






Av, iV H / 


<?cs 


At / 


Fluxes - 


= alternative 


/ kpc 


/ hours 


/ M@ 


1 Mq 


10 21 cm~ 2 


(band, 


days 


data 


name 


(ref) 


(ref) 


(ref) 


(ref) 


(refs) 


ref) 




references 


(I) 


(II) 


(HI) 


(IV) 


(V) 


(VI) 


(VII) 


(VIII) 


(IX) 



M31 r2-70 

GRO J0422+32 
= V518 Per 

LMC X-3 

A0620-00 
= V616 Mon 



784+30 
(1) 

2.49±0.30 
(2) 

50+10 

(3,4) 

1.2+0.4 

(5) 



+ 290 
120 



192 
(12) 

5.09 
(5) 

~40.8 
(13) 

7.75 
(5) 



3.97+0.95 

(15) 

~9-10 
(16) 

11. Oil. 9 

(17) 



0.46+0.31 

(15) 

~4-8 
(16) 

0.74+0.13 
(17) 



1.0+0.3, 
1.8+0.5 (12) 

0.74+0.09, 
1.6±0.4 (2,23) 



2.0 



0.76_ 

(R, 2)" 



+0.04 



0.19±0.03, 

(24-26) 



0.32 1 



1.17+0.08, 



0.58 1 



2.1 



+1.1 



(5,27) 



1+0.25 

-0.22 

(V, 33) 



1.0 



12 



0.5 2,41,42 



1.0 43,44 



33,45 



X1E J1118+480 


1. 1 1+0.05 


4.08 


6.8+0.4 


0.28+0.05 


u.uoo_ 016 , 


n kk+0.15 
u.oo_ 2g 


1.0 


6,41—43, 


= KV UMa 


(6) 


(5) 


(15) 


(15) 


0.11+0.04 (6) 


(R, 34) 




46-49 


GRS 1124-68 


5.5+1.0 


10.4 


6 0+ 1 ' 5 


0.80+0.11 


0.9+0.1 (5,28), 


0.55+0.05 




28,50-52 


= GU Mus 


(5) 


(5) 


(18) 


(18) 


-, cro + 0.42 
l- 5S -0.58 


&(B-V, 35) 






GS 1354-64 


>27 (7) 


61.1 


>7.83±0.50 


>1. 02+0.06 


2.60+0.31, 




1.0 


43,53 


= BW Cir 


(~33±6) 


(7) 


(15) 


(15) 


37.2+^ 4 (7,29) 








4U 1543-47 


7.5+0.5 


26.8 


9.4+1.0 


2.45+0.15 


1.55+0.15, 


0.68«; o i 


1.0 


43,54,55 


= IL Lup 


(5) 


(5) 


(15) 


(15) 


4.3+0.2 (5) 


(R, 36) 






XTE J1550-564 


5.3+2.3 


37.0 


10.6+1.0 


1.30+0.43 


2.5+0.6, 


0.7+0.1 


1.0 


19,56,57 


= V381 Nor 


(5) 


(5) 


(19) 


(19) 


8.7+2.1 (5) 


(V, 19) 






GRO J1655-40 


3.2+0.2 


62.9 


7.02+0.22 


2.35+0.14 


3.7+0.3 (5), 


~1.0 


1.0 


20,43, 


= Nova Sco 1994 


(5) 


(5) 


(20) 


(20) 


6.66+0.57 (20) 


(B-K, 37) 




58-60 


GX 339-4 


O+7.0 
8 -1.0 


42.1 


~5.8 


~0.52 


3.9+0.5, 


<0.3 


1.0 


43,60-64 


= V821 Ara 


(8) 


(5) 


(21) 


(21) 


61"?;? (5,30) 


(B-K, 38) 






GRO J1719-24 


2.4+0.4 


14.7 


~4.9 


~1.6 


2.8+0.6, 




0.5 


41 


= Nova Oph 1993 


(9) 


(14) 


(9) 


(9) 


41"°;°, (9,31,32) 








XTE J1720-318 










6.9+0.1, 




1.0 


43,65 


= INTEGRAL1 51 


(10) 








12.4+0.2 (10*) 








XTE J1859+226 


6.3+1.7 


9.17 


5-12 


0.68-1.12 


1.80+0.37, 


0.59+0.04 


1.0 


22,39,43, 


= V406 Vul 


(5) 


(5) 


(22) 


(22) 


8+2 (5) 


(R, 39) 




66,67 


GRS 1915+105 


9.0+3.0 


816 


14.0+4.4 


0.81+0.53 


19.6+1.7, 




1.0 


68 


= V1487 Aql 


(11) 


(5) 


(15) 


(15) 


35+3 (11) 








V404 Cyg 


4 +2 '° 


155.28 


10.0+2.0 


0.65+0.25 


3.65+0.35 


0.87+0.03 


0.5 


69-71 


= GS 2023+338 


(5) 


(5) 


(15) 


(15) 


6.98+0.76 (5,27) 


(R, 40) 







Columns give: (I) source names; (II) distance estimate; (III) orbital period of the system; (IV) mass of the compact object (black hole 
here) in solar units; (V) mass of the companion star in solar units; (VI) interstellar reddening in V-band, and interstellar HI absorption 
column (*Ay is estimated here from the relation JVh = 1-79 X 10 21 cm _2 J 4v; iPredehl fc Schmitj[T995T) : (VII) the companion star OIR 
luminosity contribution in quiescence; (VIII) The maximum time separation, At, between the OIR and X-ray observations defined as 
quasi-simultaneous; (IX) References for the quasi-simultaneous OIR and X-ray fluxes collected. References: see caption of Table 3. 



each data set, representing the average total systematic la 
errors associated with each luminosity data point. 

In addition, we searched the literature for quasi- 
simultaneous OIR and X-ray fluxes from a number of NSXBs 
in a hard X-ray state (mostly atoll sources in the 'island 
state', defined by the X-ray colour-colour diagram; se e e.g. 
lHasinger fc van der Klisll989l:lMigliari fc Fenderl2006l) . and 
HMXBs in both the hard and soft X-ray state. The same 
methodology was adopted in calculating the intrinsic OIR 
and X-ray luminosities of the NSXBs and the HMXBs. Ta- 
bles 2 and 3 list the properties of the NSXBs and HMXBs, 
respectively. A literature search for OIR Spectral Energy 



Distributions (SEDs) of BHXBs was also conducted in order 
to shed further light on the nature of the emission. Fluxes 
were used when two or more OIR wavebands were quasi- 
simultaneous. No X-ray fluxes were required for the SEDs, 
however the X-ray state of the source on the date was noted. 
Where the companion star significantly contributes to the 
emission, the estimated wavelength-dependent contribution 
of the companion was subtracted (adopting the quoted con- 
tributions given in the papers from which the data was ac- 
quired) . 
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Table 2. Properties and data collected for the 8 NSXBs. 



Source 

= alternative 

name 

(I) 


Distance 
/ kpc 
(ref) 

(II) 


Period 
/ hours 
(ret) 
(III) 


1 M Q 
(ret J 
(IV) 


M cs 
1 M Q 
(vet) 
(V) 


A v , iV H / 
10 21 cm- 2 
(refs) 
(VI) 


IJcs 

(band, 

1 CI J 

(VII) 


At 1 
days 

(VIII) 


Fluxes - 
data 

references 
(IX) 


IbK J00291+5934 


4 -0 
(72) 


2.457 
(77) 


1 . 4 

(72) 


0.039-0.160 

(72) 


1.56±0.22 
2.8±0.4 (77*) 




1.0 


97,98 


4U 0614+09 
= V1055 Ori 


3 -°-2.5 

(73) 


0.25-0.33 
(78) 


1.4 

(80) 


<1.9 (80) 
(-1.45) 


1.41±0.17 (78) 
2.99±0.01 (87) 




0.5 


99,100 


CXOU 132619.7 
-472910.8 


5.0±0.5 
(74) 








0.34±0.03 (88) 
0.9±0.1 (74) 


0.8±0.1 

(B, 88) 




74,88 


Cen X-4 
= V822 Cen 


1.2±0.3 
(75) 


15.1 
(79) 


1.3±0.8 
(15) 


0.31±0.27 

(85) 


0.31±0.16 (89) 
0.55±0.16 (90) 


0.75±0.05 (R), 
0.70±0.05 (V, 95) 


0.5 


79,95, 
101,102 


4U 1608-52 
= QX Nor 


3.3±0.5 

(5) 


-12.9 

(5) 


-1.4 
(81) 


-0.32 
(81) 


^• DO -0.18 

15±5 (14,81,91) 




1.0 


43,81 


Sco X-l 
= V818 Sco 


2.8±0.3 

(5) 


18.9 
(5) 


1.4 

(82) 


0.42 

(82) 


0.70±0.23 
1.25±0.41 (14t) 




1.0 


43,103 


SAX J1808. 4-3658 
= XTE J1808-369 


2.5±0.1 

(76) 


2.0 

(5) 


>1.7 (83) 
(-1.7) 


0.05-0.10 

(86) 


0.68±°;^ (86) 
0.11±0.03 (92) 


o o +0 - 3 

u - u -o.o 
(V, 96) 


1.0 


86,104, 
105 


Aql X-l 

= V1333 Aql 


5.15±0.75 

(5) 


19.0 

(5) 


-1.4 
(84) 


-0.6 

(84) 


1.55±0.31 (93) 
±.0±l 8 2 (94) 




0.5 


43, 

106-112 



Columns give: (I) source names; (II) distance estimate; (III) orbital period of the system; (IV) mass of the compact object (neutron 
star here) in solar units; (V) mass of the companion star in solar units; (VI) interstellar reddeni ng in V-band, a nd interstellar HI 
absorption column (*A\r and ' TVjj are estimated here from the relation Nn = 1.79 X 10 21 cm~ 2 A^: |Predehl fc Schmittlll995l) : (VII) the 

companion star OIR luminosity contribution in quiescence; (VIII) The maximum time separation, At, between the OIR and X-ray 
observations denned as quasi-simultaneous; (IX) References for the quasi-simultaneous OIR and X-ray fluxes collected. References: see 

caption of Table 3. 



2.1 Results 

Quasi-simultaneous OIR and X-ray luminosities are plotted 
of 15 BHXBs in the hard state (Fig. la), 9 BHXBs in the soft 
state (Fig. 16), 8 NSXBs in the hard state (Fig. 2a) and 9 
HMXBs (Fig. 26). W e have classed LMC X-3 as a low-mass 
XB (BHXB) because lBrocksopp et ail fcOOTal) found from 6 
years of observations on this source that its optical emission 
is dominated by long-term variations rather than the bright 
companion star in the system, due to its persistent nature 
(and we are using data from their paper). 

For BHXBs in the hard state, a strong Loir-Lx corre- 
lation exists over 8 orders of magnitude in X-ray luminosity. 
The data from all 15 individual sources lie close to this cor- 
relation but deviations in the slopes and normalisations of 
individual sources are present. The slope of the global cor- 
relation is /3=0.61±0.02, where Loir oc L^ (we do not take 
into account the Loir and Lx error bars in calculating 0). 
To calculate correlations we use the package gnuplot and ap- 
ply equal weighting to each data point. The relations will be 
biased towards sources with more data points, but we argue 
that all data points are equally important because we have 
used a maximum of one data point per day per waveband 
per source, and the sources tend to vary on timescales > 
days. 

In the soft state, the optical data (_BVJ?/-bands) lie 
close to the hard state correlation for BHXBs (Fig. 16), and 
the NIR data (JH i\-bands) lie below the correlation. A cor- 
relation is also found for the NSXBs (Fig. 2a) over 7 orders of 



magnitude in X-ray luminosity. Its slope, Loir oc /J^ 63±0 - 04 5 
is remarkably similar to the OIR-X-ray correlation of the 
BHXBs, but with a lower normalisation. Taking the two 
correlations, we find that at a given X-ray luminosity (2-10 
keV) , the OIR luminosity of a BHXB in the hard state is on 
average 19.4 times larger (1.29 dex) than that of a NSXB in 
the hard state in the X-ray range of the data. 



If we neglect the data from sources in quies- 
cence (Lx < 10 33,5 erg s _1 ), the hard state fits are 



Loir = 10 1 
10' 



13.2±0.8^0.60±0.02 



for BHXBs and Loir = 



9.7±1.8 r0.66±0.05 



for NSXBs. These are very similar to 
those with quiescent data included. The similarity of the 
fits strengthens the case for the quiescent state being an ex- 
tension to the hard state. For the analysis in the following 
Sections, we use the fits with quiescent data included, as 
they are similar enough either way. In addition, it is also 
interesting that the fit to the BHXB quiescent data alone 
has a slope Loir oc / / ^ 67±014 (there are only 4 quiescent 
NSXB data points). 

From Fig. 26 it is clear that the OIR luminosity of 
HMXBs is typically orders of magnitude larger than that 
of LMXBs, and does not appear to correlate with Lx- This 
is consistent with the high-mass companion star dominating 
the OIR emission. The large range in OIR luminosities be- 
tween sources is likely to be due to the differing masses and 
spectral types of the companions. Since this interpretation 
agrees with the evidence in the literature, we will not discuss 
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Table 3. Properties and data collected for the 9 HMXBs. 



Source = alternative name 
(I) 


Compact 
obj Get 

(ii) 


Distance 
/ kpc (ref) 
(HI) 


A v (ref) 
(IV) 


N n 1 10 21 cm' 2 

("ref") 

(V) 


At 

j days 
(VI) 


Fluxes - 

data references 

(VII) 


SMC X-3 


NS 


58 1+5 6 H 1 31 


1 5+n 7 (120) 


2 9+1 4 ("120"! 




130 131 


CI Cam = XTE J0421+560 


unknown 


5±2 (114) 


3.2+1.2 (121) 


5+2 (124) 


1.0 


43,124,132 


LMC X-4 


NS 


50+10 (3,4) 


0.31+0.06 (122*) 


0.55+0.10 (122) 




43,133 


A0535+26 = HDE 245770 


NS 


2t° A 7 (115) 


2.3+0.5 (115) 


11.8+1.5 (125) 


1.0 


125 


GX 301-2 = 4U 1223-62 


NS 


5.3+0.1 (116) 


5.9+0.6 (116) 


20+10 (126) 




116,126 


V4641 Sgr = SAX 1819.3-2525 


BH 


9.6+2.4 (5) 


1.0+0.3 (5) 


2.3+0.1 (127) 


0.2 


43,134,135 


KS 1947+300 = GRO J1948+32 


NS 


10+2 (117) 


3.38+0.16 (117) 


34+30 (117) 


1.0 


43,117 


Cyg X-l = HD 226868 


BH 


2.1+0.1 (118) 


2.95+0.21 (123) 


6.21+0.22 (128) 


1.0 


43,136-138 


Cyg X-3 = V1521 Cyg 


unknown 


10.3+2.3 (119) 


20+5 (119) 


85+1 (129) 


1.0 


43, 139 



Columns give: (I) source names; (II) BH = black hole, NS = neutron star; (III) distance estima te; (IV) interstellar reddening in V-band 
(*Ay is estimated here from the relation iVjj = 1.79 X 10 21 cr?i~ 2 Av; [Predehl fc Schmitt 1995); (V) interstellar HI absorption column; 
(VI) The maximum time separation, At, between the OIR and X-ray observations denned as q uasi-s imultaneo us; (VII) Ref erences for 
the quasi-simultaneous OI R and X-ray fluxes colle cted. Re f erence s for Tables 1, 2 fc 3: (1) IStanek fc Garnavichl .19981): ( 21 
iGelino fc Harrison! .2003lj : OIlBovd et al.| .20001) : ( 41 iKovacj .20001) : (5llJonker fc Nelemansj 120041) : (6llChatv et al.H2003l) : (7) 



ICasares et alj [20041): ( 8) IZdziarski et a" 
IChapuis fc CorbeJ [20041) : (12l lWmiams" 




Kolbl [2003ft; 

Hvnes et alJ 



(251 IWilms et alj 



) | Cowlev et al 
' (211 lHvnes et al 




91 Idella Valle. Mirabel fe Rodriguez! il994ft : (101 ICadolle Bel et alj||200 



et al 




; (13^^^^^^^^^^^ 
, (17) Gclino, Harrison 

... ._ [2003aft : C2 2llHvnes*et al.l <2002al) : (23llShrader et all [1997ft : (241 ISoria et al, 

; (26) IWang et alj .20051): (271 iKong et alj .20021): (28llEbisawa et al] .19941): (291 iKitamoto et al 



14l | Liu. van Paradiis fcwmd 
[200lft; (18)lEsin et al] [200Cft 



den Heuvell 12001,. 
( 191 lOrosz et al] 



Jjrr, 



(30) Zdziarski et al. (1998); 



es et alj .20041): f35llOrosz et al] .1999) 
. (39llZurita et all 12002ft ; (40 1 ICasares et al 



(311 iTanakal [19931: (32 
(36llOrosz et al 



(441 [ Brocksopp et al 

[200Cft: (49 1' 
ISut aria et al] I 
_J2C 



(33) McCli ntock Horne : , 



R.emillardMl995l): (34~ 



T TOrosz et al] .19981): (37llGreene et al] .200 j): (381 IShahbaz. Fender fc Charles! 
1.1993) : (411 brocksopp et al] [2004ft : (42llGarcia et al] J200lf) ; (4 3) RXTE ASM ; 



; (45 ) iKuulkerj .19981): (461 iMcClintock et al. | .2003): (471 iKiziloglu et al] .200 j): (481 lllemura et alj 



12005) 



(2003): 



(57) .1 
Fender! 



.20021): (62 


iHoman et al] {200511: (631 


J20031: (671 


lHaswell et all 1200011: (68) | 



[Corbel 
iTomsick et al 

.19991); (711|Han fc Hiellmind .1992D: (7 2 
IChevalier et alJ ll98Sll: (76l|in't Zand et alj 1200 

^8^^^^^^^^^^j_et_^uE^^jr7^1^__^_^_^_™_-^- 

: Yound [200d): (85llTorres et a* 



ml 



(501 iKing. Harrison fc McNamaral ]l99g): (51 1 Idella Valle 
Brocksopp ct al .1 .2001bl): (54l | Buxton fc Bailvnl 



20Mi5- 
llTm 



._, Masctti & Bianchim 

iKalemci et al] <200Rl): (561 



199811: (521 



Hame urv et al] 
" l2002j); ( 61) 



2005j^( 59l lTorres et al] .20051): (60l lChatv et al] _ 

[20041: (64l llsrael et al] .20041): ( 65l lNa.gata et al] .2003): (66 1 
l2000D : (69llHvnes et alj|2004h: f70llZvcki. Done fc Smith! 



Gallow ay et alj .2 0051: (73) Bran dt et al 
(771 IShaw et al] -200511 



(1992); (74) Rutlcdgc ct al. (200 



_; ( 821 ISteeghs fc Casaresl ._ _ 

(86llWang et al] ]200ll) : (8 71 ISchulJ .1999ft 



752 

(78llNelemans et alj 120041): (791 ICampana et all 12004 
" 12002ft ; (8 3llCampana et alJ <2002}i~?84) 

lHaggard et al]' 12004 



IWelsh^ 

(89llBlair et al] ]l984l): (901 iR.utledge et alJ l200lf): (911 iGrindlav fc Lillerl ]l97^i 7l ( 921 ICampana et all 20051) : (931 IChevalier etal] 
<1999ft: (941 ICampana fc StellaH2003ft : (951 IShahbaz. Navlor fc Charles! <1993ft : (961 iBurderi et alj .2003ft: (971 mares et al] l200i " 



(981 ISteeghs et al.1 120041) : (991 iMachin et alJ ]l990l) ; ( 1001 Russell et al. (in pr eparati on); (1011 iKaluzienski. Holt fc~ 
(102llCanizares. McClintock fc Grindlavl ]l980ft: 



198! 



(106l lCharles et al J 11981 
( lllllMaitra et alj 12004 



McNamara et al] .2003ft: (10 41 ICampana fc Stellal .2004ft : ( 105 1 iHomer et all 
iOj); (107) iJain et al] ]l999ft: (1081 [Jain et al] l200Cft: ( 109l lMaitra et alj i2003ft: (HOllMaitra et af] 
Mbft ; (ll^pMaitra fc BailvrJ J20031: (113llColel [l998li: (114llMioduszewski fc Rupenl i2004): ( 1151 



1998); (116) Kapcr et al] ]l995ft: (1171 INcgucrucla et al] .2003): (1181 iMassev. Johnson fc Dcgioia-Eastwood 



Ivan Kerkwiik et al] .19961): (1 201 iLequeux et al] 1199211: (12 11 iHvnes et al] .2002bft: (122llNaik fc Paul [20031: (1231 [Wi 

(124llOrlandini et al] ]200(i): (1251 lOrlandini et al.1 pool): (126' 
ISchulz et al] [2002ft: (1291 IStark fc Saial [20031): (1301 IClark et alj 



Mukhcricc 



(131 



PauJll2004l): (1271 iDickev fc Lockmar 
Cowley fc Schmidtkel [20041) : (132llClark n 



(1331 iHeemskerk fc van Paradiis! .198^1: (13 41lK ato et al.1 jl 9*9Sft1 (1351 iBuxton. Bailvn fc _Maitral[2005rJl7 (1361 
f): (1371 iBrocksopp et alj 



1999rJh : (1381 IBrocksopp et al] [l999aft : (1391 iKoch-Miramond et"al 




Bochkarcv . 
' [2002ft 



(119) 
11982); 
128) 
2000ft: 



Lvutvi 



the HMXBs in further detail (although they are mentioned 
in Section 3.5). 

There are clear advantages and disadvantages of these 
OIR-X-ra y compilations com pared to the I/ ra dio _ ix ap- 
proach of iGallo et ail l|2003ft . The hard state OIR-X-ray 
correlation in BHXBs includes data from many BHXBs in 
low-luminosity quiescent states (7 sources with 10 30 < L x < 
10 33 5 erg s" 1 ), which was not possible for the radio-X-ray 
corre lation due to radio detector limits (but see lGallo et alJ 
2006). The BHXB sample also includes sources in the LMC 



and M31, which are too distant (and hence faint) to observe 
at radio wavelengths. However, unlike the radio-X-ray com- 
parison, HMXBs cannot be included in these OIR-X-ray 
plots as the OIR emission is dominated by the companion. 
In addition, two sources situated in the galactic plane (IE 
1740.7-2942 and GRS 1758-258) were included in the radio- 
X-ray correlation but are not observable at OIR wavelengths 
due to the high levels of extinction towards the sources. 

In Section 3 we attempt to interpret the relations found 
between OIR and X-ray luminosities in terms of the most 
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Figure 3. Spectral Energy Distributions (SEDs) of 15 BHXBs. The key in the top left panel corresponds to the X-ray state of the source 
on the date of observation. A key to the slope of the continuum, i.e. the spectral index a (where L v oc u a ), is indicated in the second 
panel on the top row. The dates of observations and references for the data are listed in Table 4. 
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Table 4. Dates and references for the data in Fig. 3. 



Source 


X-ray 


Dates X-ray 


Dates 


X-ray 


Dates 


References 




state 


(MJD) state 


(MJD) 


state 


(MJD) 






hard 


48866, 48870, 49006,7, 49039, 49209, 49590 






1-4 


IMP Y_5 
i_J 1V1 K_j y\ O 


hard 


50151, 50683 soft 


50324, 51038 


unknown 




5-7 


A0620-00 


hard 


42792, 43097-100, 42859 soft 


42650-2, 42703 






8 - 11 


XTE J1118+480 


hard 


51652, 51975, 53385-6, 53393, 53409, 53412 






12, 13 2 


GRS 1124-68 


hard 


48453, 48622 soft 


48367-71 






14 - 16 


GS 1354-64 


hard 


50778, 50782, 50851, 50888 








17 


4U 1543-47 


hard 


52486, 52490, 52495, 52501 soft 


52454, 52469 






18 


XTE J1550-564 


hard 


51260, 51630, 51652, 51717 soft 


51210, 51660, 51682 






19, 20 


GRO J 1655-40 


hard 


53422 soft 


50217, 50254, 50286, 50648 




21 - 23 


GX 339-4 


hard 


44748-59, 49539, 50648, 52368, 52382 




intermediate 


52406 


23 - 25 


GRO J1719-24 


hard 


49096, 49101, 49108, 49266, 49280, 49539 








23, 26 - 29 


XTE J1720-318 


hard 


52782 soft 


52659, 52685, 52713 






30 


XTE J1859+226 


hard 


51605-8 soft 


51465,9, 51478, 51488, 51501 




31, 32 


GS 2000+25 1 


soft 


47360 unknown 


47482 






33, 34 


V404 Cyg 


hard 


47684, 47718-27 soft 


47678 






35 - 38 



1 This so urce (alternative name: QZ Vul) is not tabulated in Table 1; its distance and reddening are 2.7±0.7 kpc and Av ~ 3.5, 
respectively ( Jonker & Nclemans 2004). 2 This previously unpublished data was obtained with the Liverpool Telescope and the United 
Ki ngdom Inf rared Telescope. See Br ocks opp et al. (in preparation) for th e data redu ction recipe used. Refere nces: (11 [Bartolini et al] 
" 19941 ): ( 21 iGoranskii et al.Ul99rt: (31 ICastro-Tirado. Ortiz fc Gallegol <1997fr: (4llHvnes fc: HaswelH <199flh: (51 iBrocksopp et alJ 




l2001aD: (6 1 ITreves et al] ll987lW71 ITreves et al. <198Sft : ( 81 [Robertson. Warren fc Bvwaterl [l97 
IOkelil977h: I'll) iKleinmann Brecher fc Inghaml Il97fitl : ( 121 Ichatv 
(14) iKing et al] <1996ft: (15l ldella Valle et al] Il998tl : ( 16llBailvnl ]l 



(i97n^^^t^ ^ 1 <2001all: (20llJain et al.l fcOOlbll: 

l2002il: (241ICorbel fc Fenderl <2002ft: (25llHoman et al] 



„ Greensteinl il97' 
et alT l200^1~7l3l this paper (see B rocks opp et al, 
1992)1: (171 IBrocksopp et al.U2001rJ): (18llBuxton , 



(10) 

in preparation ); 
: Bailvr] <2004h; 



(211 iBuxton Bailvn fc Maitral i2005al1: (2 2 1 iHvnes et al] Jl998ft: (231 IChaty et al] 
— i2005l): (261 ISekiguchi fc van Wvkl <1993tl: (271lAllen fe Jones! Il993j: (281 
lAllen fe Gilmorei ]l993ft: (291 IBrocksopp et al] <2004l); (30llNagata et al] <2003l); (311lHaswell etUd] l2OO0tl: (32llHvnes et al] <20o"2af 
(33l lCharles et al] Il98lh : (341 IChevalier fc Ilovaiskvl 1 199011; (351 ISzkodv et al] il989T>: (361 iGehrz. Johnson fc Harrisonl <1989l) : (371 

ICasares et al.Ml99ll) : (381 lHan fc Hiellmind jl9 92li 



likely dominant emission processes. OIR SEDs were col- 
lected from 15 BHXBs in a range of luminosities and X-ray 
states and are presented in Fig. 3 and the dates of obser- 
vations and references are given in Table 4. The SEDs are 
interpreted in Section 3.3. In Sections 3.4 and 3.5 we dis- 
cuss additional patterns, applications and implications of 
the empirical relations. The results and interpretations are 
summarised in Section 4. 



3 INTERPRETATION & DISCUSSION 

3.1 Jet Suppression in the Soft State 

The significant drop in some of the OIR data in the soft state 
for BHXBs compared to the hard state (Fig. 16) is not due 
to changes in the X-ray luminosity during state transition 
because although the X-ray spectrum significantly changes 
during transition, the bo lometric (and 2-10 keV) X-ray lu- 
mino sity does not (e.g. IZhang et al] Il997t iKording et al] 
2006). The jet component at radio wavelengths is known 
to decrease/increase in transition to/from the soft X-ray 
state , due to the quen ching of the jets in the soft state 
(e.g. iGallo et al]l2003t) . In Fig. 4 we split the OIR data 
from BHXBs into three wavebands: BV R-bands, /-band and 
JHK -bands, and plot the monochromatic OIR luminosity 
(L„; i.e. flux density scaled for distance) against Lx, over- 
plotted wit h the radio data L Vtla .dio-Lx for BHXBs from 
IGallo et al] (2003; 2006). We see that the normalisations in 
L v (as well as the power-law slopes) for the radio-X-ray and 
OIR-X-ray hard state correlations are similar to one another 



for BHXBs: at a given Lx, the radio and OIR monochro- 
matic luminosities are ~ equal, implying a flat spectrum 
from radio to OIR wavelengths for all BHXBs in the hard 
state. We return to the hard state interpretation in Section 
3.2. 

Fig. 4 also shows a clear suppression of all the JHK 
data, and a little, if any of the BVR in the soft state. The I- 
band data appear to sit in the centre of the two groups. All 9 
BHXBs for which we have soft state data are consistent with 
this behaviour, suggesting it is ubiquitous in BHXBs. We 
interpret this as the NIR wavebands being quenched as the 
jet is switched off, as is observed at radio wavelengths. The 
NIR appears to be quenched at a higher X-ray luminosity 
than the radio, but this effect may be just due to an upper 
Lx limit adopted bv lGallo et al.l i2003fl when compiling their 
radio-X-ray data. 

The optical wavebands in the soft state lie close to the 
OIR hard state correlation. Most of the optical data lies 
below the centre of the hard state correlation, with the ex- 
ception of V404 Cyg, whose optical luminosity is enhanced 
in the soft state. There is still a debate as to whether V404 
Cyg entered the soft state, and the data here were taken 
close to the supposed state transition. The 7-band appears to 
be the "pi vot" point, as already shown bv lCorbel fc Fenderl 
(2002) and Homan et al. (2005) to be where the continuum 
of the optically thin jet meets that of (possibly) the disc. The 
NIR quenching in the soft state implies that this waveband 
is dominated by the jets in luminous hard states, just be- 
fore/after transition to/from the soft state. The optical data 
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hard state; BVR-bands 
hard state; I-band 
hard state; JHK-bands 
hard state; OIR fit 
soft state; BVR-bands 
soft state; I-band 
soft state; JHK-bands 
both states; Radio (GHz) 



30 31 32 33 34 35 36 37 38 39 
log(L x ; ergs" 1 ; 2-10 keV) 

Figure 4. X-ray 2-10 keV luminosity versus OIR and radio monochromatic luminosities for BHXBs. 



are not quenched, suggesting a different process is dominat- 
ing at these wavelengths. 

An alternative interpretation could be that the disc 
dominates the OIR in both the hard and soft states, but 
changes temperature during state transition, shifting the 
blackbody from (e.g.) the OIR in the hard state, to the 
optical-UV in the soft state. This would have the effect of 
reducing the NIR in the soft state but maintaining the opti- 
cal, as is observed. We argue that this is not the case because 
there is evidence in many BHXBs for two spectral compo- 
nents of OIR emission, t he redder of which is quenched in the 
soft state (Section 3.3: | Jain et alJl2001bt iBuxton fc Bailvr] 
l2004lHoman et alJl2005tl . 

The jets should contribute negligible OIR flux in the 
soft state, so we can estimate the OIR contribution of the 
jets at high luminosity in the hard state from the level of soft 
state quenching. The mean offset of the OIR soft state data 
from the hard state fit is 0.30±0.32 dex, 0.71±0.21 dex and 
1.10±0.26 dex in L^om for the BVR, I and JHK-bands, 
respectively. This corresponds to a respective fractional jet 
component of 50^gQ, 81~t 13 and 92±5 percent in the BVR, 
I and JHK-bands. It is clear that the spectrum of the jet 
extends, with spectral index a ~ 0, from the radio regime 
to the NIR in BHXBs. The position of the turnover from 
optically thick to optically thin emission must be close to 
the NIR waveband for the radio-NIR spectrum to appear 
flat (unless the optically thick spectrum is highly inverted, 
which i s not seen in the radio spectrum most of the time; 
see e.g. lFendeJl200l . The OIR spectrum in the hard state 
is flatter than in the soft state, as is expected if the jet 



component is present in the former and not in the latter; 
this is explored further in Section 3.3. 



3.2 Modelling Disc and Jet Contributions in the 
Hard State 

We now attempt to interpret the empirical hard state 
correlations in terms the three most cited OIR emis- 
sion processes: X-ray reprocessing in the disc, the vis- 
cously heated disc and jet emission. We adopt the the- 
oretical model between op ti cal an d X-ray luminosities of 
Ivan Paradi is & McClintoc kl Jl994l) for X -ray reprocessing 
(Section 1.1). Ivan^aradTsj^McClintoc^ normalised their 
relation to a sample of systems, some of which are BHXBs 
and some NSXBs. Here, we normalise the relation to our 
sample of NSXBs (just the op tical -BV^-RI-bands), which is 
much larger than the sample of lvan Paradiis fc McClintocld . 
There is more concrete evidence for NSXBs to be domi- 
nated by X-ray reprocessing in the optical regime than in 
BHXBs (see Section 1). We also take into account a depen- 
dency of the relation on sys tem mass that was neglected by 
Ivan Paradiis fc McClintockL but is needed here to compare 
between NSX Bs and BHXBs. From Keple r's third law and 
equation 5 of lvan Paradiis fc McClintocld (and the discus- 
sion that follows), we have: 

Lopt oc L][ 2 a oc L][ 2 {M co + M CS ) 1/3 P 2/3 

In the left panel of Fig. 5 we plot the monochromatic 
luminosity (optical and NIR) versus L^/ 2 a for each BHXB 
and NSXB data point in the hard state (adopting the val- 
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Figure 5. BHXBs and NSXBs in the hard X- ray state, with models for X-ray r eprocessing (left) and jet (right) components. The X-ray 
reprocessing models are derived from vPfcKf_jyagi_Paradii s h. McClintoclJIf a nd the jet models assum e a flat spectrum from radio 
to OIR and use the relations from GFP lOallo et al.ll200.Tl for BHXBs and M&F iMigliari fc FenderlEoOrl) for NSXBs. 



ues of P, M co and M cs for each source from Tables 1 and 
2). The solid line represents the power-law fit to the opti- 
cal NSXB data, fixing the slope at unity (Lopt oc L^/ 2 a). 
We expect a lower normalisation for the NIR data due to 
the shape of the OIR spectrum, which may have a spectral 
index 0.5 < a < 2.0 (a conservative range based on theo- 
retical and empirical results of X-ray reprocessing; see e.g. 
lHvneslEo05h . The upper and lower dotted lines indicate the 
expected correlations of the NIR data (approximating the 
optical to the l/-band centred at 550nm and the NIR to the 
.ff-band at 1660nm) if a — 0.5 and 2.0, respectively. The 
dashed line shows the fit to the optical BHXB data (fixing 
the slope at Lopt oc L^[ 2 a), which appears to be elevated 
by ~ 1 order of magnitude in L^oir compared to the op- 
tical NSXB data. Similarly, the NIR data of both BHXBs 
and NSXBs lie above the expected correlation for emission 
from X-ray reprocessing. The implications of these fits are 
discussed below. 

For the jets, we take the models described in Sec- 
tion 1.1 where the spectrum is flat from radio to OIR, 
and norm alise them u s ing th e em pirical L Ta< g -Lx re l ations 
found by iGallo et alJ l2003t) and iMigliari fc Fender! fc006h 
for BHXBs and NSXBs, respectively. In the right panel of 
Fig. 5 we plot Z/„,oir versus Lx for hard state BHXBs and 
NSXBs. The same relation is expected between optical and 
NIR jet emission when plotting monochromatic luminosity 
because we are assuming the jet spectrum is flat; a ~ 0. 
We find that the jet models can approximately describe the 



optical and NIR data of the BHXBs. The NIR data from 
NSXBs lie above the expected jet correlation but possess a 
similar slope, whereas the slope of the optical NSXB data is 
very different to that predicted from jet emission. 

The models for OIR emission from a viscously heated 
disc are described as follows. For the simplest viscously 
heated steady-sta te disc, there are two limiting regimes 
jFrank et al.ll2002l) . For hv < kT, we expect Loir oc m 1/4 . 
This is simply the Rayleigh- Jeans limit and will only apply 
well into the IR. For hv 3> kT, the relationship is steeper, 
Loir oc rh 2//3 . For typical disc edge temperatures of 8,000- 
12,000 K, expected power-law slopes (Loir oc rh 7 ), are cal- 
culated to vary from 7 ~ 0.3 in the Tf-band, to 7 ~ 0.5 in 
V , and 7 ~ 0.6 in the UV. Using the calculations linking 
Lx and rh in Section 1.1, this corresponds to expected cor- 
relations of the form Loir oc L^, where 0.15< /3 <0.25 for 
BHXBs and 0.30< (3 <0.50 for NSXBs. 

A summary of the results of fitting these models to the 
observed data is provided in Table 5. It is clear that for 
BHXBs, the slopes j3 of the observed relations can be ex- 
plained by the X-ray reprocessing model or the jet model, 
but not by the viscous disc model (Homan et al. 2005 also 
ruled out a viscous disc origin to the OIR emission in the 
hard state of GX 339-4 because the mass accretion rate in- 
ferred from the luminosity is much higher than expected). 
This is also true for the NIR data of the NSXB sample, 
however the slope of the NSXB optical data cannot be ex- 
plained by the jet model and are accurately described by 
both the reprocessing and viscous models. At high luminosi- 
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Table 5. Parameters for the hard state models: theory versus observations. 



Sample X-ray reprocessing model Jet model — Viscous disc model — 

Model | /3data-/3model I * Model | /3 data -/3 modcl | Model | /3 data -/3 modc i | 

"-model "model 



BHs: I/j^opt 


nL^a 


0.05±0.03 


9.3±0.4 




0.11±0.02 


1.05±0.07 


rO.25 

L x 


0.34±0.02 


BHs: Z/„,nir 


( VNIR ) a nL° 5 a 


0.06±0.03 


15.5-81. 3 1 " 




0.09±0.04 


1.78±0.16 


r0.17 
L x 


0.44±0.04 


NSs: L„,opt 




0.09±0.02 


1.0* 


rl.4 
L X 


0.81±0.03 


6.03±1.94 




0.09±0.03 


NSs: L u , Nm 


( VNIR ) a nL°- 5 a 

u O PT 


0.05±0.03 


3.2-16.6 1 " 




0.09±0.41 


9.55±3.08 




1.19±0.41 



For | Pdata-/3modei |, /3 and n are free parameters; for ™ da ' d a l , P is fixed at the value of the model and n is a free parameter; * 
n is defined by the fit to the optical NSXB data (see text); ' the range corresponds to an OIR spectral index 0.5 < a < 2.0. 



ties in the NSXB sample, no OIR data points lie far below 
the jet relation, suggesting that this process may indeed play 
a significant OIR ro le here (as is seen in a few NSXBs, see 
Migli arT et al .11200(1 and references therein). 

The normalisation n of the BHXB data is closer to the 
jet model than the reprocessing model. Although this is con- 
sistent with the constraints derived from Section 3.1 (the 
jets are contributing ~ 90 percent of the NIR luminosity 
here), an optical excess of ~1 order of magnitude from the 
reprocessing model is not expected. The excess is unlikely 
to be fully explained by the jets because of the lack of op- 
tical quenching in the soft state (Section 3.1). In addition, 
the optical spectrum of most BHXBs is inconsistent with jet 
emission dominating (Section 3.3). Instead we suggest that 
OIR emission from reprocessing is enhanced for BHXBs at 
a given Lx due to the localisation of the source of X-rays. 
For example, the X-ray emitting region in BHXBs may have 
a larger scale height than in NSX Bs and will therefore illu - 
minate the disc more readily fe.g. iMiniutti fc Fab ian 2004). 
This may account for the high value of n for the optical 
BHXB data, but still struggles to explain the even higher n 
for the NIR BHXB data. 

We note that there are deviations of individual sources 
from the correlations which may be caused by distance and 
reddening errors, or by differing emission process contri- 
butions due to the range of orbital separations or slight 
differences in the slope of the radio-OIR jet spectrum 
between sources (or other system parameters not consid- 
ered). For example, XTE J1118+480 has a small disc and 
is known to produce significant optical jet emiss ion (e.g. 
iMalzac. Merloni fc Fabianl2004t |Hvnes et aill200d) whereas 
V404 Cyg possesses a large di sc and is dominated by X- 
ray reprocessing in the disc (e.g. lWaener et allll99ll) . Other 
OIR emission mechanisms that could contribute, for exam- 
ple disc OIR emission due to magnetic reconnection, have 
not been modelled here and could also contribute to the 
scatter in Fig. 5. These processes cannot be ruled out, but 
are unlikely to easily explain the observed correlations. 

In addition, the orbital inclination may affect the level 
of OIR emission, in particular from X-ray reprocessing 
in the disc. To explore this, we have plotted in Fig. 6 
the average normalisation n of the data for each source 
against the known estimated orbital inclination, fixing the 
slope at Lorn oc L^ 6 for each source. The best estimate 
of the inclina tion of each s ystem (where known) wa s 
obtained from iRitter fc Kolbl (l200St iTreves et alJ lll986t 
iFomalont. Geldzahler^^BTadshawi ll20rjlF ! Iwan^^t^l] 
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Figure 6. Mean normalisation n of each source, fixing the slope 
at Loir 0(1 ^x 6 ' versus orbital inclination. 

<200ll) and iFalanga et alJ J20051) . We find no evidence for 
there being a direct relation between Loir and the orbital 
inclination i, suggesting the effect of inclination is subtle. 

3.3 The Spectral Energy Distributions 

Although no clear patterns are visible from the SEDs in Fig. 
3 on first inspection, closer analysis reveals support for many 
of the conclusions made so far. Optically thin synchrotron 
emission is expected to produce an OIR spectral index a < 
0, and this is seen in part of the SEDs of 10 out of 14 BHXBs 
in the hard state. In the soft state only 3 out of 10 BHXBs 
are observed to have a < 0, suggesting a synchrotron com- 
ponent playing a larger role in the hard state than in the soft 
(a in some sources could be dominated by uncertainties in 
A v ; e.g. the SEDs of XTE J1550-564 are red even in the soft 
sta te, which is likely due to an underestimated extinction; 
see iJonker fc Nelemanslliooi) . In comparison, a > is seen 
in the SEDs of 10/14 BHXBs in the hard state and 10/10 
in the soft state. These spectra are likely to have a thermal 
origin and agree with re cent analysis of optical/UV SEDs of 
6 BHXBs iHvnesll200t) . The SEDs generally appear redder 
in the hard state than in the soft, as is observed by the NIR 
suppression in the soft state (Section 3.1). 

A clear suppression of the NIR, and not the optical 
bands in the soft state is visible in th e SEDs of GX 339- 
4, XTE J1550-564 and 4U 1543-47 llHoman et alJ l2005l : 
I Jain et alll2001bt iBuxton fc Bailvnll2004ll . In contrast, the 
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value of a in XTE J1720-318 and XTE J1859+226 ap- 
pear not to change between the hard and soft states. The 
SEDs show no substantial evidence for the turnover in the 
jet spectrum from optically thick (a ~ 0) to optically thin 
(a ~ —0.6) synchrotron emission, as we may expect to see 
in th e hard state (except in GX 339-4; see lCorbel fc Fendeil 
2002). This is consistent with the turnover lying close to 
the NIR (Section 3.1), possibly just redward of K-band. In 
some systems, a is more negative at low luminosities. This 
effect could be the result of (a) a cooler accretion disc, or (b) 
the jets contributing more than the disc at low luminosities. 
Since the former process cannot explain the steeply negative 
SEDs of a few BHXBs at low luminosities, we suspect both 
processes may play a role. 

Finally, the hard s t ate NIR excess see n in GX 339-4 
JCorbel fc Fended Eool iHoman et alJl2005l : interpreted as 
where the optically thin jet spectrum meets the blue thermal 
spectrum) is here also seen in LMC X-3 (tentatively), 4U 
1543-47, XTE J1550-564 and V404 Cyg. In these sources, 
we interpret the NIR excess as originating in the optically 
thin part of the jet spectrum. It is interesting to note that 
OIR SEDs of a number of NSXBs show an IR excess that 
cannot be explained by thermal emission , and is likely to 
originate in the jets iMigliari et al ] |2006t and references 
therein). 

3.4 The Rise and Fall of an Outburst 

During a transient outburst typical of BHXBs, the source 
will either remain in the hard state for the entire out- 
burst (e.g. XTE J11 18+480; GRO J1719-24; see also 
iBrocksopp et all l2004Tl or transit into the soft (or in- 
termediate or very high) state, before returning to the 
hard state and declining in luminosity. A hysteresis ef- 
fect has been identified in BHXB transients, and at 
least one NSXB transient: the hard to soft state tran- 
sition always occurs at a higher X- ray luminosity than 
the t r ansition back to the hard state l^ accaronefc Coppil 
20031: iMaccarone Gallo fc Fended 120031: lYu et alJ l2003f 
Fender et alJl2004l : IHoman fc Bellonill2005lb 
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Here, we present a prediction of the OIR behaviour of an 
outburst that results from this hysteresis effect and the mod- 
els in Section 3.2. The following sequence of events should 
occur for a BHXB outburst that enters the soft state: 

• Lorn and Lx increase in the hard state rise of the out- 
burst, with jets and reprocessing contributing to the OIR. 

• The source enters the soft state, quenching the jet com- 
ponent. Loir drops but Lx is maintained. 

• Loir and Lx decrease before transition back to the hard 
state, when the jet component returns. 

• Loir and Lx continue to decline, with jets and X-ray 
reprocessing contributing towards Loir. 

This sequence is illustrated in the top panel of Fig. 7. 
In this schematic, we fix the hard-to-soft state transition 
at L x ~ 10 38 erg s" 1 and the soft-to-hard, at L x ~ 10 37 
erg s . To test our hysteresis prediction, we split the hard 
state BHXB data into data from the rise of outbursts and 
data from the declines. We define rise and decline as before 
and after the peak X-ray (2-10 keV) luminosity of the out- 
burst, respectively. We do not include data on the decline 
for sources that remained in the hard state throughout the 
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Figure 7. Top panel: Schematic of the expected Z/oiR^^x be- 
haviour for a BHXB outburst that enters the soft state. Lower 
panel: As the top panel but with OIR data from the rise, soft 
state and decline (for sources that entered the soft state) stages 
of outbursts. 



outburst, as this is not what our prediction is testing. In 
the lower panel of Fig. 7 we plot these and the soft state 
data. Data from persistent sources (LMC X-3 and GRS 
1915+105) and from sources in or near quiescence are not 
included. We find that the prediction is consistent with the 
data, with some inevitable scatter from errors as described 
in Section 3.2. The loop is larger in the NIR data, as is ex- 
pected because the jets contribute more in the NIR than in 
the optical regime. Other reasons for any deviations from 
the expected models are also discussed in Section 3.2. 



3.5 Applications of the Correlations 

The existence of the OIR-X-ray correlations leads to a 
number of intriguing tools and uses for quasi-simultaneous 
multi-wavelength data. 

The Mass Accretion Rate: 

In Section 1.1 we show how Lx and rh are thought to 
be linked for BHXBs and NSXBs. Here, we can use the 
empirical hard state Loir-Lx correlations to link the OIR 
luminosity to the mass accretion rate: 

BHXBs: Loir oc L£ 6 oc m 12 
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NSXBs: Loir oc L x ' 6 oc m 06 

Essentially, Loir for BHXBs and NSXBs respond to 
Lx in the same way, but not to rh, since Lx varies with rh 
differ ently for BHXBs and NSXBs. From equations (1) and 
(7) of Kor ding et alJ i2006T) we can estimate accretion rates 
directly from Loir in the hard state: 

BHXBs: Lom/erg s' 1 « 5.3 x 10 13 (m/g s" 1 ) 1 ' 2 

Or: m/g s" 1 « 3.7 x 10" 12 (Lom/erg s" 1 ) 8 

NSXBs: Lom/erg s" 1 « 3.2 x 10 23 (m/g s -1 ) ' 8 

Or: m/g s" 1 « 6.7 x 10- 40 (Lom/er fl s" 1 ) 17 

Given the level of the scatter in the correlations, we 
expect these calculations to be accurate to ~ one order 
of magnitude. The OIR luminosity is all that is required 
to estimate rh for hard state objects, however rh would 
be more accurately measured from Lx, where most of 
the energy is usually released. In addition, it is possible 
to estimate the X-ray, OIR and radio luminosities in the 
hard state quasi-simultaneously, given the value of just one 
because they are all linked through correlations. 

Parameters of an X-ray Binary: 

Quasi-simultaneous OIR and X-ray luminosities can con- 
strain the nature of the compact object (BHXB or NSXB), 
the mass of the companion (HMXB or LMXB) and the 
distance and reddening to an X-ray binary. This is possible 
because the data from BHXBs (in the hard and soft states), 
NSXBs and HMXBs lie in different areas of the Lx-Loir 
diagram, with some areas of overlap. If the distance and 
reddening towards a source is known (not necessarily at 
a high level of accuracy), its quasi-simultaneous OIR and 
X-ray fluxes could reveal the source to be any one of the 
above types of XB. In addition, if the nature of the compact 
object is known (BH or NS), but the distance and/or 
reddening is not, the fluxes can constrain these parameters. 
We stress that the total errors associated with the data 
(top left error bars in each panel of Figs 1 and 2) need to 
be considered to define the areas of overlap. 

Current techniques used to infer the nature of the com- 
pact object in XBs include X-ray timing analysis (e.g. ther- 
monuclear instabilities on accreting neutron stars produce 
Type I X-Ray bursts), the X-ray spectrum (the well-known 
X-ray states of BHXBs and tracks in the colour-colou r di- 
agrams of NSXBs; e.g. iMcClintock fc Remillardll2006tl and 
optical timing analysis in quiescence (the orbital period and 
radial velocity amplitude constrain the mass function) . This 
new tool has the power to constrain the nature of the com- 
pact object requiring only Loir, Lx and, only at high lumi- 
nosities, the X-ray state of the source at the time of obser- 
vations. 

This tool may have many applications. X-ray all-sky 
monitors such as the RXTE ASM are continuously discov- 
ering new XBs which are subsequently identified at optical 
wavelengths. In addition, campaigns are underway to find 
cxtragalactic XBs, many of which have optical counte rparts 
discovered with the HST (e.g. IWilliams et alJ I2005T) . It 
would also be interesting to see where ULXs and SMBHs 
(e.g. the X-ray and NIR flares seen from Sagittarius A*; 
e.g. lYuan. Quataert fc NaravarJ 12004) lie in the Loir-Lx 
diagram with the inclusion of a mass term, and whether this 



could be used to constrain the BH mass in these systems. 
The Power of the Jets: 

The evidence in this paper confirms that in the literature of 
the presence of a flat jet spectrum from radio to NIR wave- 
lengths in the hard state, with the optically thin-optically 
thick turnover close to the LT-band (see Sections 3.1 and 
3.3). The power of the jets is sensitive to the position of the 
turnover since it is dominated by the higher energy photons. 
However, the mo s t recent calculations of the jet power (e.g. 
i Gallo et alJl2005l : iHeinz fc Grimmll2005t IMigliari fc Fended 
2006) already assume the jet spectrum extends to the IR so 
no extra calculations are required from this work. 



4 CONCLUSIONS 

We have collected a wealth of OIR and X-ray fluxes from 
33 XBs in order to identify the mechanisms responsible 
for the OIR emission. A strong correlation between quasi- 
simultaneous OIR and X-ray luminosities has been discov- 
ered for BHXBs in the hard state: Loir oc L x - 61±0 ' 02 . This 
correlation holds over 8 orders of magnitude in Lx and in- 
cludes data from BHXBs in quiescence and at lar ge distances 
(LMC and M31), which were (until recently; see lGallo et alJ 
2006) unattainable for radio-X-ray correlations in XBs. All 
the NIR (and some of the optical) BHXB luminosities are 
suppressed in the soft state; a behaviour indicative of syn- 
chrotron emission from the jets at high luminosities in the 
hard state. A similar correlation is found for NSXBs in the 
hard state: Loir oc / / ^ 63±0 04 ] w hich holds over 7 orders of 
magnitude in Lx. At a given X-ray luminosity, a NSXB is 
typically 20 times fainter in OIR than a BHXB. 

Comparing t he hard state OIR data to radio data of 
IGallo et alJ <l2003l) . we find that the radio-OIR jet spectrum 
in BHXBs is ~ flat (F v = constant) at a given Lx (see 
also lFendeill200 ill . From this and OIR SEDs of 15 BHXBs, 
which show an IR excess in a number of systems, we deduce 
that the optically thick-optically thin turnover in the jet 
spectrum is likely to be close to the A"-band for BHXBs. In 
comparison, th e turnover pro bably lies further into the IR 
for NSXBs f see IMigliari et aUl2006ft . 

By comparing the observed OIR-X-ray relations with 
those expected from models of a number of emission pro- 
cesses, we are able to constrain the mean OIR contributions 
of these processes for XBs. Table 6 summarises the results. 
We find from the level of soft state quenching in BHXBs that 
the jets are contributing ~90 percent of the NIR emission at 
high luminosities in the hard state. The optical BHXB data 
could have a jet contribution between zero and 76 percent 
but the optical SEDs show a thermal spectrum indicating 
X-ray reprocessing in the disc dominates in this regime. In 
BHXBs, ambiguity arises from the fact that the slope of the 
expected OIR-X-ray relations from the jets and X-ray re- 
processing are essentially indistinguishable. Emission from 
the viscously heated disc may contribute at low luminosi- 
ties in BHXBs, but cannot account for the observed cor- 
relations. In the NSXBs, the correlations can be explained 
by the X-ray reprocessing model alone (agreeing with the 
current thinking), but the jets may play a role at high lu- 
minosities, especially in the NIR. The exact contributions 
of the emission processes are likely to be sensitive to many 
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Table 6. The OIR emission processes that can describe the em- 
pirical OIR-X-ray relations and SEDs. 



Sample X-ray X-ray Jet Viscous Intrinsic 

state repro- emis- disc com- 
cessing sion panion 



NSXBs; OPT 


hard 


V 


X 


V 


X 


NSXBs; NIR 


hard 


V 


V 


X 


X 


BHXBs; OPT 


hard 


V 


V 


X 


X 


BHXBs; NIR 


hard 


X 


V 


X 


X 


BHXBs; OIR 


soft 


V 


X 


V 


X 


HMXBs; OIR 


all 


X 


X 


X 


V 



individual parameters, such as the size of the accretion disc 
and the shape of the jet spectrum (e.g. the OIR luminosity 
of the jets is very sensitive to the slope of the optically thick 
jet spectrum). 

The SEDs show a non-thermal component (a < 0) in 
most sources in the hard state, and thermal emission (a > 0) 
is present, probably in all sources in both hard and soft 
states. The soft state could originate from a combination of 
X-ray reprocessing ( as the OIR-X-ray re lations suggest) and 
the viscous disc (e.g. lHoman et al.l2005|) . The SEDs of many 
BHXBs are redder at low luminosities, which seems to be 
due to both a cooler disc blackbody and a higher fractional 
contribution from the jets. BHXBs may be jet-dominated 
at low luminosities (for more on jet-dominated states, see 
| FendeJ200H:lFalcke. Kording fe Markofi1l2004l iFender et"ai] 
120041: iGallo et alJl2005t iKording et al.ll200d) . We have also 
made a prediction of the 'LoiR-ix path' of a BHXB for a 
typical outburst, based on a hysteresis effect whereby the 
hard state rise reaches a higher luminosity than the hard 
state decline, if the source enters the soft state. The data 
obtained in this paper appear to agree with the prediction, 
with some inevitable scatter. 

Since the X-ray, OIR and radio luminosities and the 
mass accretion rate are all linked through correlations in the 
hard state, it is possible to estimate the quasi-simultaneous 
values of all of these parameters, given the value of just one, 
by e.g. daily monitoring of the X-ray or OIR fluxes (which 
is currently being done for some sources by the RXTE ASM 
and by ground-based telescopes). In addition, we have dis- 
covered a potentially powerful tool to complement current 
techniques, that can constrain the nature of the compact ob- 
ject, the mass of the companion and the distance/reddening 
towards an XB, given only the quasi-simultaneous X-ray and 
OIR luminosities. Data from BHXBs (in the hard and soft 
states) , NSXBs and HMXBs lie in different areas of the -Lx- 
Lom. diagram, with small areas of overlap. The tool is most 
useful for e.g. faint sources with poor timing analysis or at 
large distances, i.e. extragalactic XBs and new transients. 

Further work that could constrain the emission pro- 
cess contributions include identifying linear polarimetry 
perpendicular to the jet axis (a diagnostic of optically 
thin synchrotron emission), emission line equivalent width 
analysis (in particular testi ng for the Baldwin effect; e.g. 
iMushotzkv fc Ferlandlll984h and analysing Loir-Ioc rela- 
tions in the intermediate and very high X-ray states and in 
ULXs and SMBHs. 
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